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ABSTRACT 

In the last decades we witnessed an increase in studies of open clusters of the Galaxy, especially because of the good determination for 
a wide range of values of parameters such as age, distance, reddening, and proper motion. The reliable determination of the parameters 
strongly depends on the photometry available and especially on the U filter, which is used to obtain the color excess E(B-V) through 
the color-color diagram (U-B) by (B-V) by fitting a zero age main-sequence. Owing to the difficulty of performing photometry in the 
U band, many authors have tried to obtain E(B-V) without the filter. But because of the near linearity of the color-color diagrams that 
use the other bands, combined with the fact that most fitting procedures are highly subjective (many done "by eye") the reliability 
of those results has always been questioned. Our group has recently developed, a tool that perfoms isochrone fitting in open-cluster 
photometric data with a global optimization algorithm, which removes the need to visually perform the fits and thus removes most 
of the related subjectivity. Here we apply our method to a set of synthetic clusters and two observed open clusters (Trumpler 1 and 
Melotte 105) using only photometry for the BVRI bands. Our results show that, considering the cluster structural variance caused only 
by photometric and Poisson sampling errors, our method is able to recover the synthetic cluster parameters with errors of less than 
10% for a wide range of ages, distances, and reddening, which clearly demonstrates its potential. The results obtained for Trumpler 1 
and Melotte 105 also agree well with previous literature values. 
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1. Introduction 

Our research group has a particular interest in using open clus- 
ters in our studies of the spiral structure of our Galaxy. Some 
interesting results were obtained with data of these objects 
that are compiled in our catalog (| Dias et al.(2002)) : the rotation 
speed of the spiral pattern as 25 ± 1 km^\s^'kpc^^ and the co- 
rotation radius at 1 .06 + 0.08 x the solar Galactocentric distance 



(Dias & Lepine(2005) I and Lepine et al.(2008) determined the 



' epicyclic frequency to be 4-3 + 5 km 's ^kpc ' at the solar ra- 
dius. 

In this kind of study and in many others that involve open 
\ clusters it is crucial to accurately determine the fundamen- 
. tal parameters (distance, reddening, and age) of the objects. 
Traditionally this is done using the color-magnitude diagrams 
via the main-sequence fitting (MSF) technique, which in the ma- 
■ jority of cases is fit to the data by visual inspection. The work 
' of |Paunzen & Netopil (2006 ) (hereafter PN06) indicates that the 
usual method appUed up to now can determine distances of open 
clusters with absolute errors of less than 20% (see all details in 
the cited paper). 

From the observational point of view, the determination of 
fundamental parameters of open clusters requires accurate pho- 
tometry (errors lower than 0.03 in magnitude and lower than .05 
in color index), deep photometry (mainly for objects not previ- 
ously studied) and also seems to require data in the U-filter The 
U-filter in principle allows a better determination of the color 
excess and therefore leads to an accurate determination of the 
distance modulus via MSF. The determination of the color ex- 
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cess toward the cluster is performed based on the observed pho- 
tometric data, typically by visual main-sequence fitting of the 
(B-V) vs. (U-B) diagram. Results on determining the color ex- 
cess without the use of the U-filter are unreliable if done in the 
usual manner by visual fits because of the near linearity of the 
color-color diagrams and the subjectivity of the analysis. 

We have addressed the subjectivity of the isochrone fits in 
Monteir o et al.(2010^ (hereafter MDCIO), where we developed 
a global optimization tool using the cross-entropy optimization 
algorithm (CE) to automatically fit open cluster UBV photomet- 
ric data with well defined criteria (see paper for details). We 
showed that the method was robust and the results obtained for 
10 open clusters agreed well with previous studies found in the 
literature. Although the work used only UBV photometry and 
the usual sequence to fit isochrones (determine E(B-V) and then 
distance and age), the results obtained with the global optimiza- 
tion indicated that it might be possible to address the color excess 
determination without the use of the U-filter The advantages of 
such a procedure are clear, given that the U-filter is possibly one 
of the most problematic ones to get good data for, which typi- 
cally requires longer exposure times. 

In this work we investigate the possibility of determining 
open cluster fundamental parameters from BVRI photometry 
alone. In Sec. 2 we briefly describe the global optimization 
method developed in MDCIO and the adaptations made for the 
present problem. In Sec. 3 we describe benchmark tests using a 
set of synthetic open clusters generated with known parameters. 
In Sec. 4 we apply the method to the well studied open clusters 
Trumpler 1 and and Melotte 105 and discuss the results. In Sec. 
5 we present our conclusions. 
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2. The cross-entropy optimization method 



3. E(B-V):fromO.Oto3.0. 



The CE procedure provides a simple adaptive way of estimat- 
ing the optimal model parameters. Basically, the CE method in- 
volves an iterative procedure where each iteration consists of: 

1 . Random generation of the initial parameter sample, respect- 
ing pre-defined criteria; 

2. Selection of the best candidates based on some mathematical 
criteria; 

3. Random generation of updated parameter samples from the 
previous best candidates to be evaluated in the next iteration; 

4. Optimization process repeats steps (ii) and (iii) until a pre- 
specified stopping criterion is fulfilled. 

In MDCIO we fully describe the algorithm and the fitting 
method applied to open clusters and we refer the reader to that 
work and references therein for more details. In short, the algo- 
rithm generates a set of possible isochrone solutions (simulated 
open clusters) given a pre-defined initial mass function, binary 
fraction, and metallicity. Each generated solution is then com- 
pared to the observed data through an objective function. The 
best solutions are selected, a new set of solutions is generated 
from them and the process iterated until convergence. In the 
present work we maintain the complete structure of the previ- 
ous algorithm discussed in detail in MDCIO, changing only the 
objective function to use the BVRI data instead of only UBV. 

In the case of isochrone fitting, the objective function used is 
the likelihood of the data (in the present case BVRI photometry) 
for a given model isochrone, which is obtained from the product 
of probabilities of the m stars observed and is given by: 



P{BVRI\ISOmod)^Yj 



nO-BCTvO-RO-j 



j^^p_ {B-B,so) ^ 



EXP- 



Icte 
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To simplify the analysis we kept the metallicity constant at 
the solar value. This should have no major impact on the evalua- 
tion of the capability of the method to determine the fundamental 
parameters because the changes due to metallicity are relatively 
small when compared to the other parameters. 

3. Synthetic clusters 



1.00 - 

0.95 - 

0.90 L_ 
1.10 

1.05 7 

1.00 - 

0.95 - 



Errors Related Age 



E(B-V) 



log(age) 



log(age) 

Fig. 1. Relative errors for fits of the synthetic clusters fitted with 
our method and BVRI photometry only as a function of synthetic 
cluster age. The synthetic clusters have fixed E(B - V) - 0.4, 

jcr phot 

= 0.8% and are at a distance of 1000 pc. 



where B, V, R, and I are the observed magnitudes and cr^, 
cry, ctr, and ctj are the photometric uncertainties in each filter 
Here ISOmod designates the model isochrone for a set of pa- 
rameters (E(B-V), distance and age) and Bjso, Vjso, Riso^ and 
I ISO are the model magnitudes for each point generated from the 
model. 



X = ]~[p(ByW5<9„„^), 



(2) 



where X is the likelihood obtained in the usual manner. 

As in MDCIO, the tabulated isochrones are taken from 
[Girardi et al. (2000)| and [Marigo et al (2008) | and are specified 
by four parameters, namely, age, distance, extinction constant 
and metallicity. In the isochrone fitting we defined the parameter 
space as follows: 

1. Age: from logiage) - 6.60 to log(age) - 10.15; 

2. distance: from 1 to 10000 parsecs; 



To test the capability of recovering fundamental parameters 
using the CE method we constructed a series of synthetic open 
clusters. The generated clusters were chosen to be representative 
of clusters present in the catalog ( |Dias et al.(2002)) and w ere 
generated from the isochrone tables of Girardi et al. (2000)] and 
Marigo et al (2008)1 which are the same as those used in the fit- 
ting algorithm. 

We also generated clusters with distinct numbers of stars 
from low to high density to evaluate the capacity of the method 
to retrieve the parameters in unfavorable conditions. 

All clusters were generated adopting solar metallicities, sam- 
pled from a Salpeter IMF and with a 50% binary population at 
a distance of 1000 pc, E(B - V) = 0.40, 3o-pho, = 0.8% and 
limiting my magnitude of 18. All fits were executed 20 times to 
perform an error analisis via bootstrap re-sampling. 

We did not include simulated contamination from field stars, 
but this evaluation is already under way. But as we show in 
MDCIO, the filtering method we developed removed a good part 
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Errors Related to Number of Stars 



Errors Related to E(B-V) 
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Fig. 2. Relative errors for fits of the synthetic clusters fitted with 
our method and BVRI photometry only as a function of number 
of generated stars. The synthetic clusters have fixed E{B -V)- 
0.4, 3o-pi,gi = 0.8%, log(age) - 8.0 and are at a distance of 
1000 pc. 



Fig. 3. Relative errors for fits of the synthetic clusters fitted with 
our method and BVRI photometry only, as a function of number 
of generated E{B - V) values. The synthetic clusters have fixed 
logiage) - 8.70 at a distance of 2100 pc. 



of the contamination and we expect that the main effect would 
likely be an increase in the estimated uncertainty for the param- 
eters. 

The results obtained for the synthetic clusters fitted with our 
method and BVRI photometry only are summarized in Fig. 1 and 
Fig. 2 where the relative errors for fits as a function of synthetic 
cluster age and relative errors for fits as a function of number 
of generated stars are shown. The distance and E(B - V) were 
not varied in the results presented in Figs. 1 and 2 because the 
major changes were caused by the variation in age and number 
of stars. Varying the distance only affects the photometric errors, 
which we assumed here to be some percentage of the observed 
magnitudes, and so would only increase the final uncertainty. In 
Fig. 2 the age was also fixed at logiage) - 8.0. 

To check the reliability of the method under distinct red- 
dening values, we performed a set of tests with synthetic clus- 
ters generated with E{B - V) varying from 0.3 to 3.0, covering 
the observed parameter range, and with logiage) - 8.70 yr at 
a distance of 2100 pc. We note that an increase in the magni- 
tude limit to mv - 22 was necessary for the synthetic clusters 
with EiB - V) > 2 so that the main sequence could be mini- 
mally sampled. The results presented in Fig. 3 indicate that the 
uncertainty of the distance determination tends to increase with 
increasing EiB - V), which is to be expected. Interestingly, the 
uncertainties obtained for the various reddening values studied 
were essentially constant around 0.03, which shows up as de- 
clining relative errors in Fig. 2. 



We also tested the fitting method on a synthetic cluster gen- 
erated with the same parameters as the one presented in Sec. 5.1 
of MDCIO, i.e., 400 stars with photometric error 'icTphoi - \% 
sampled from an isochrone with logiage) - 8.70 yr, at a distance 
of 2100 pc, EiB - y)=0.40 and Z=0.019. The results obtained 
from the fitting algorithm were 

- tog(flge) = (8.73 ± 0.05) yr 

- distance — (2100 ± 70) pc 

- EiB - y) = (0.39 ± 0.03) 

The results indicate that the method using only BVRI is ca- 
pable of recovering the cluster parameters with about the same 
accuracy as the UBV method presented in MDCIO. The preci- 
sion is also about the same exept for the reddening, where the 
B VRI-fitted value shows larger uncertainty. 

4. Application to observed clusters 

We applied this new approach of fitting isochrones to a cou- 
ple of well-studied open clusters. The selected clusters were 
Trumpler 1 and Melotte 105, which have observational data in 
the BVRI filters obtain ed by |Yad av fe Sagar (2002) (hereafter 
YS02) and Sagar, Mun ari, & de B oer (2001) (hereafter SMBOl) 
respectively. The main reason for choosing these particular clus- 
ters is that they are the only ones that have BVRI data and were 
also studied in our previous work (MDCIO), where we applied 
the CE method using UBV photometry only, following the tradi- 
tional route of first determining the reddening through the (U-B) 
versus (B-V) diagram and then fitting the theoretical isochrones. 
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Table 1. CE fit and literature results for Trumpler 1 



Ref 



E(B - V) Distance (pc) log(Age(yr)) 



PN06 0.57 : 

MDCIO 0.59 : 

YR02 0.60 : 

this work 0.60 : 



0.04 
0.05 
0.05 
0.03 



2356: 

2419: 
2600: 

2706: 



511 
185 
100 
: 120 



7.5 ±0.1 
7.8 ± 0.2 

7.6 ± 0.1 
7.2 ± 0.3 



Table 2. CE fit and literature results for Melotte 105 



Ref 



E(B - V) Distance (pc) log(Age(yr)) 



PN06 0.48 ± 0.05 

MDCIO 0.47 ±0.02 

SMBOl 0.52 ± 0.05 

this work 0.42 ± 0.03 



2094 ± 159 
1750 ± 111 
2300 ± 200 
1932 ± 95 



8.4 + 0.1 
8.4 ±0.1 

8.4 ± 0.1 

8.5 ±0.1 



The results obtained for Trumpler 1 and Melotte 105 using 
our fitting method and the BVRI photometry agrees very well 
with the results obtained by YS02 and SMBOl and the results 
from our previous work using only UBV photometry. The val- 
ues also agree with those of PN06. In Fig. 4 we show the final- 
fit results (isochrones) for Trumpler 1 and in Fig. 5 for Melotte 
105, overplotted on color-magnitude diagrams where the sym- 
bol sizes are proportional to the weight determined for each star 
based on the filtering scheme described in MDCIO. 

For Trumpler 1 the results obtained from our BVRI fit agrees 
well with the values obtained by YS02 (see Table 1). Some small 
but significant differences are found when we consider the re- 
sults of PN06 and MDCIO, but these are likely due to weights 
given to the bright stars as discussed in MDCIO. The weight of 
the bright stars are Ukely the main cause for the age differences. 
Notice, however, that using the BVRI only we obtain almost the 
same E(B-V) as in the fit made with the (U-B) versus (B-V) di- 
agram. 

For Melotte 105 we find a similar situation where we found 
small significant differences in E(B-V) and distance when com- 
paring our results with those of SMBOl. The differences are 
again likely due to distinct weights given to particular stars. The 
weight might be particularly relevant when comparing the re- 
sults of MDCIO, where the red giants received a lower weight 
compared to this work, which could be the reason for the lower 
distance encountered in that work. 



5. Conclusions 

We investigated the possibiUty of fitting theoretical tabulated 
isochrones to open cluster BVRI photometry with a global op- 
timization method developed in MDCIO. In the context of up- 
coming ground-based telescopes and space missions that will 
produce massive amounts of data, this tool, wich additionally 
has the possibility of full automation, is of great value. 

The results obtained with synthetic open clusters that are 
generated with parameters that cover typical values found in the 
Uterature show that even without the use of the U-filter, funda- 
mental parameters can be determined with a relative precision 
of about 10% or better in most cases in the studied parameter 
space, also with an accuracy similar to the method described in 
MDCIO using UBV photometry only. In particular we showed 
that the reddening can be determined for a wide range of clus- 
ter ages, despite the lack of visual structure in the color-color 
diagrams of the clusters. Even when considering objetcs with 
few stars the method was able to recover the original parameters 



within the uncertainties, showing some systematic behavior for 
clusters with a low number of stars as shown in Fig. 2. 

It is important to point out that the relative uncertainties dis- 
cussed previously are to be taken as lower limits to the real ex- 
pected uncertainties that are attainable with real data. The er- 
rors reported are only due to the adopted error in the photometry 
and Poisson statistics that arise from drawing a finite number of 
stars from model isochrones. Real clusters will have structural 
variations because of binary fraction or chemical composition 
among many other factors, will affect the final uncertainty of 
the method. Because it is difficult to quantify the nature of these 
factors and their effect on the final uncertainty, a definitive error 
analysis should be carried out in a case-by-case basis. 

Applying the method to two observed clusters also stud- 
ied with our optimization method and UBV photometry (see 
MDCIO for details) showed that the use of BVRI filters only also 
produced consistent results. The results agree well with those 
obtained by previous works, in particular when considering the 
results of PN06, where the authors included a large number of 
independent literature determinations. Assuming that the results 
of PN06 are the best estimates of the real cluster parameters, we 
can estimate a better value of the final uncertainty of the method 
for these objetcs. For Trumpler 1 the highest relative uncertainty 
comes from the distance and is about 15%. For MeUote 105 the 
largest relative error is due to the reddening being approximately 
13%. These values are likely more realistic estimates of the real 
uncertainty of the method for these two clusters. 

Our results show that it is indeed possible to reliably deter- 
mine open-cluster fundamental parameters from BVRI photom- 
etry only despite the lack of visual features in the color-color 
diagrams. 
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Fig. 4. Final-fit results for Trumpler 1 overplotted on color-magnitude diagrams where the symbol sizes are proportional to the 
weight determined for each star as described in MDCIO. 
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Fig. 5. Same as Fig. 3 for Melotte 105. 



